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CHAPTER  1 
INTRODUCTION 


1. 1 Background 

Characteristic  modes  have  long  been  used  In  the  analysis  of  radiation 
and  scattering  by  dielectric  and/or  magnetic  bodies  whose  surfaces  coincide 
with  coordinate  surfaces  in  coordinate  systems  for  which  the  Helmholtz  equa- 
tion is  separable.  From  consideration  of  the  scattering  matrix,  Garbacz  [1] 
has  shown  that  similar  modes  must  exist  for  any  material  body.  An  extensive 
theory  for  perfectly  conducting  bodies  was  given  in  reference  [1],  but  the 
dielectric  and  magnetic  body  case  was  not  developed.  An  alternative  treat- 
ment of  the  characteristic  modes  for  perfectly  conducting  bodies,  starting 
from  the  impedance  operator  for  the  conducting  surface,  has  been  given  by 
Harrington  and  Mautz  [2].  The  computation  of  such  modes  has  also  been  con- 
sidered by  Harrington  and  Mautz  [3].  A theory  of  characteristic  modes  for 
dielectric  bodies,  magnetic  bodies,  and  for  bodies  both  dielectric  and 
magnetic,  has  been  developed  by  Harrington,  Mautz,  and  Chang  [4].  In  this 
work,  a theory  of  characteristic  modes  for  material  bodies  is  developed 
using  equivalent  surface  currents.  This  is  in  contrast  to  the  approach  used 
in  [4],  which  used  the  induced  volume  currents. 

The  modes  are  defined  by  a weighted  eigenvalue  equation  in  such  a way 
that  both  the  generalized  network  matrix  [5]  and  the  scattering  matrix  [1], 

[2]  for  the  body  are  diagonalized.  The  presentation  given  in  this  work  leads  to 
explicit  formulas  for  determining  the  mode  currents  and  fields  of  two-dimensional 
objectB.  The  formulas  remain  the  same  for  dielectric  bodies,  magnetic  bodies,  and 
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for  bodies  both  dielectric  and  magnetic.  In  particular,  the  scattering 

problem  of  a two  dimensional  material  cylinder  will  be  presented.  This 

formulation  of  the  problem  is  applicable  to  any  general  material  body. 
Details  are  worked  out  only  for  two-dimensional  problems. 

1-2  The  Fundamental  Operator  Equation 

Let  the  material  body  be  represented  as  in  Figure  l-l. 


Figure  1-1.  A general  material  body. 


jjf , *>  incident  fields  (Wavy  underline  denotes  vector 

- inside  fields  quantities.) 

E , H - outside  fields 


The  problem  of  Figure  1-1  can  be  viewed  as  a linear  superposition  of 
two  cases, 


(1)  zero  field  inside 
(II)  zero  field  outside. 


These  two  cases  are  illustrated  in  Figure  1-2. 


n 
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Figure  1-2.  A decomposition  of  the  original  problem. 


In  caae  (I),  let  eQbe  the  material  constants  in  the  zero-field 
regiont  and  similarly  in  case  (II) , let  u*  e be  the  material  constants  in 
the  zero-field  region.  Having  done  so,  radiation  formulas  [6]  for  unbounded 
space  can  be  employed.  The  J and  M are  equivalent  surface  currents  [6].  Since 
there  are  no  actual  surface  currents,  the  following  conditions  should  be 
satisfied  by  the  equivalent  currents 


£ 


M'  + M" 

A*A 


0 

0 


(1-1) 

(1-2) 


Equations  (1-1)  and  (1-2)  come  from  the  fact  that  tangential  components  of 
fields  are  continuous  across  the  interface  in  the  original  problem.  Note 
that 

E - - jwA*  - 71'  - - V x F1  + E1  (1-3) 

-O  AAA  6 E AA*  AAA  A/‘SA 

W Aaa  O 

H - - Ju>5"  -71'+-7  xA'  +H1  (1-4) 

AAA  m A*"  AAA  AAA 


and 
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(1-5) 
(1-6) 

(1-7) 
(1-8) 

vector  potential  due  to  electric  current 

F - vector  potential  due  to  magnetic  current 

scalar  potential  due  to  electric  charge 

d - scalar  potential  due  to  magnetic  charge 
m 

Primed  quantities  refer  to  case  (1) 
doubly  primed  quantities  refer  to  case  (II) 

In  terms  of  general  operator  notations,  the  following  equations  are  obtained 


E - - jwA"  - V - ~ 7 * F" 

VA  <*VS  ^ AAA 

H - - juF"  - V$"  + - 7 x A" 

W AA.?1  AAA. 

and 

E8  - - jwA'  - 7$'  - i-  7 x F' 

AAA  AAA  aaa?  E0  AAA.  AAA 

H8  - - jwF'  -7*'  + -7*A' 
J in  p 

//A  ^ 0 AAA  AAA 

I 8 

where  E And  H are  scattered  fields,  and 

AAa  /s#A 
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where  the  definitions  of  the  operators  are  obvious  when  comparing  equations 
(1-9)  and  (1-10)  with  (1-5)  to  (1-8).  Note  that  the  tangential  fields  are 
continuous  at  the  boundary  surface,  l.e. 


n x (E8  + E1  - E)  ■ 0 

/S/vS  /*-\  aA*.  AAA 

d-11) 

n x (HS  + H1  - H)  - 0 

/*«A  AAA  A/A  AAvS 

where  n - unit  normal  pointing  outward,  or 

(1-12) 

8 1 
- n x E + nxE**nxE 

AAA  /VA  A/VS  A'VN  /*\\  V* 

(1-13) 

-nxH8  + n-<H«*nxH^ 

AM  AM  A* A AAA  AAA 

(1-14) 

The  following  equations  are  obtained  by  substituting  equations  (1-9)  and 


(1-10)  into  (1-13)  and  (1-14). 

- n x l'(J')  - n x L'  (M’)  + n * L''(J")  + n x L''  (M")  - nxE1  (1-15) 

- n X L'  (Jf)  - n X Li,(M')  + n x I"  (J")  4-  n x l"  (M")  - nxH1  (1-16) 

KM  a/sm  ms  MV.  /vs*  “nw  A/A  AAA 

By  equations  (1-1)  and  (1-2),  it  follows  that 

- n x [(L'  + L"  )J']  - n x [(Li,  + L"  )M' ] - n x El  (1-17) 

/WV  11  11  Ass  m/sm  11  12  AA  ^ 

- n X [<L!  + L")J']  - n x [(L'  + L"  )M’ ] - n x H1  (1-18) 

AA\  iL  AV,  LL  ,yv>  A/SA  AAS 

Define  all  the  operators  to  be  tangential  operators;  the  above  two 


equations  can  be  put  into  standard  matrix  form  as: 


(1-19) 


■<lh  + Lll)  -<Li2  + Lll> 


<l21  + L21> 


■<L22  + L22> 


m m 

J' 

E1 

Mf 

tan 

H1 

AM 

AM 

tan 


The  [ ] means  the  tangential  components  of  the  bracketed  quantity  on  the 
boundary  surface.  Let 


Le'  ■ - Le"  ■ 


„ tM  t • « „ 7 * T •*  m — T ** 

Lll*  m L22«  Lm  L22 


and  define 


- C' (M’ > - - — V x p«  - L*  (M') 

'VV'V 


- C"(M')  - - 7 x F"  - L"  (M') 


C'(J’)  - *7*  V x A'  - U.  (J’) 

/WS 


C'CJ')  - “ 7 x A"  ■ I,"  (J') 

A*A/*V  ^ am  A/'.N 


Hence  equation  (1-19)  becomes 


Le'  + Le" 


■(C'  + C") 


C'  + C" 


Lm'  + Lm" 


J' 


M' 

tan 


H 


tan 


(1-20) 


It  Is  convenient  to  rearrange  equation  (1-20)  Into  the  form 
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► i 

I f 


where 


Le  -jC  J 

-jC  Ln  JM 


Le  - Le'  + Le"  J - J' 
Lm  ■ Lm'  + Lm"  M « M' 
C - C'  + c" 


(1-21) 


and  the  subscript  "tan"  has  been  dropped  for  brevity.  Equation  (1-21)  is 
simply  the  familiar  operator  equation  expressed  below. 


L(f)  - g 

L Is  a tangential  operator  on  the  surface  of  the  material  body,  and 


(1-22) 


Le  -jC 
-jC  Lm 


5 g “ 


1-3  Format 

In  this  work  the  impressed  magnetic  field  is  assumed  to  be  axially- 
directed  (perpendicular  polarization).  The  derivation  of  all  the  formulas 
for  an  impressed  axially-directed  electric  field  (parallel  polarisation) 
will  not  be  giver,  in  the  main  body  of  this  work,  ln  order  to  conserve  space, 
however,  explicit  formulas  will  be  provided  in  Appendix  A.  A list  of  computer 
programs  will  be  given  in  Appendix  B. 

The  content  of  this  work  1b  as  follows.  In  Chapter  2,  the  operator 
equation  is  reduced  to  a matrix  equation  suitable  for  numerical  computation. 
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The  reduction  is  accomplished  by  using  the  method  of  moments  for  perpendicular 
polarization.  The  equivalent  surface  currents  can  be  obtained  by  matrix  inver- 
sion. A concise  theory  of  characteristic  modes  for  material  bodies,  based  on 
the  surface  formulation,  and  explicit  formulas  for  obtaining  the  modal  solutions 
are  given  in  Chapter  3.  Chapter  4 is  a presentation  of  calculate  t made  for 
cylinders  of  different  material  constants  using  both  the  matrix  inversion  method 
and  the  modal  method.  Chapter  5 la  a discussion  of  the  results.  The  computa- 
tions presented  in  this  work  were  performed  on  an  IBM  System  370,  Model  155 
digital  computer.  The  computer  programs  are  written  in  FORTRAN  IV. 
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CHAPTER  2 

MATRIX  FORMULATION 


The  determination  of  equivalent  surface  currents  requires  the  solution 
of  the  following  Inhomogeneous  equation 


L(f)  - g (2-1) 

where  L is  the  matrix  of  operators 


and 


Le 

-JC 

•jc 

Lm 

f - 


r 

J 

J 

g " 

Ju 

4 

JM 

aaa 

JH 

- 

(2-2) 


(2-3) 


This  chapter  presents  the  reduction  of  equation  (2-1)  to  matrix  form  by  the 
method  of  moments. 


2. 1 Method  of  Moments 

To  apply  the  method  of  moments , an  appropriate  inner  product  for  the 
problem  is  ( indicates  transpose) 

<f,  g>  ■ If  f g ds 

AA*  J ) *** 

8 

- (J  • E - H • M)ds  (2-4) 

t ' A//,  A/V>  W' 

8 

A solution  by  the  method  of  moments  is  obtained  as  follows.  Define  electric 
expansion  and  testing  functions  as 
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I 


e n 
m A*A 


and  magnetic  expansion  and  testing  functions  as 


The  expansion  for  f la  then  of  the  form 


f - T (if'  + vf®) 

L n n n n 

A">  n MSA  (M* 

where  the  I and  V are  constants  to  be  determined. 
n n e 

The  Inner  product  of  each  {w  } with  equation  (2-1)  yields 


<W*  L(f  )>  - <W*  g> 
m n m . ^ 


where 


<w\  L(f  )>  - <W*  l (I  L(f*)  + V L(f“)> 
m n m * n n nn 

'***  ***  ^ n vA  MSA 


2 I <M*.  L(f®)> 
_ n ta  n 

n AW  AAA 


+ 2 V <H*,  l(f“)> 
n in  n 

n msk 


• Le(J  )ds  + 
n 


><;•  *>  ■ fj  *c  • e1j* 

■*a  “A  " AAA  AAA 


(2-10) 
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Similarly, 


and 


Equation 


vhera 


the  inner  product  of  each  {W^}  with  equation  (2-1)  yields 

<W“,  L(f  )>  - <W“,  g>  (2-11) 

in  n in 

aw  aaa 


<w“,  L(f  )>  - 
in  n 
w A** 

+ 

<<->  - {{ 

1-8)  and  Equation  (2-11)  can  be  placed  in  matrix  form 


(Z] 

[B] 

[I] 

[v1] 

[C] 

1* 

[V 

tv] 

m 

[i1] 

I Tn  1J  K ■ «««.>*■ 

n | AlV.  /VV*> 

I v JJ  w“  • L (M  )de 

“ n ' ' m m n 
n g 

w“  • j^ds 


(2-12) 

(2-13) 


(2-14) 


Z - 

mn 


an 


an 


Y - 

an 


t r 


Le(J  )ds 
n 

AAA 

(-jC) (Mn)de 

AAA 

<-jC)(J  )ds 
n 

AAA 


L (M  )d« 
m n 


(2-15) 


(2-16) 


(2-17) 


(2-18) 
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V1- 

m 


We  * E^ds 
m 

am  ma 


(2-19) 


l1  - 
m 


<C  • J»ld* 


(2-20) 


Choose  J ■ W®.  Note  that  [Z]  is  obviously  symmetric,  already  shown  by 
n n 

AAA 

Herrington  and  Mauts.  With  the  choice  H ■ w , [Yj  la  the  dual  of  [Z], 

n n 


magnetic  instead  of  electric,  so  the  synmetrlc  nature  of  [Y]  can  be  easily 
established.  It  is  known  that  C(M  ) gives  rise  to  an  electric  field  and 


C(J  ) will  produce  a magnetic  field.  Observe  that  by  reciprocity 
n 


[ (E*  • Jb  - H*  • Mb)  - ff  (Eb  ♦ J*  - Hb  • M*)d 

J A/NA  AAA  aaa  J AAA  AAA  AAA  AAA 


(2-21) 


Now,  consider 

(i)  In  situation  "a"  only  electric  sources 
(ii)  In  situation  "b"  only  magnetic  sources. 
It  follows  that 


f[(-H*  • M^de  - [(  <Eb  • Ja)d« 

))  AA*  )) 


(2-22) 


Hence 


B - C 
mn  run 


Consequently,  [C]  is  the  transpose  of  [B],  or 


(2-23) 


IB]  - [C]  with  J - W* 
n n 


M - W® 
n n 


(2-24) 
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To  this  point  the  matrix  formulation  is  completely  general  and  has 
been  achieved  without  reference  to  specific  excitation,  expansion  functions, 
end  tasting  functions.  Mote  that  every  one  of  the  operators  I.e,  C,  and  Lm 
is  composed  of  two  parts  as  indicated  by  equation  (1-21) , and  consequently 
every  matrix  element  in  [Z],  [B],  [C]  and  [Y]  has  two  parts;  one  is  due  to 
the  primed  operator,  and  the  other  is  due  to  the  doubly  primed  operator. 

2.2  Expansion  Functions  and  the  Evaluation  of  [Z]  Matrix  Element 

In  this  section,  the  Incident  field  is  the  axially  directed  magnetic 
field,  H*.  Before  going  into  any  specific  excitation  for  the  scattering 

NVS 

problem,  some  general  considerations  about  the  evaluation  of  different  types 
of  matrix  elements  will  be  presented  as  follows. 

Note  that  the  original  problem,  the  scattering  from  material  bodies, 
has  been  decomposed  into  two  cases,  and  their  associated  operators  are  of  the 
same  functional  form.  F or  instance,  the  expression  for  Le"(J)  will  be  identl- 
cal  to  Le'(jJj  except  for  the  constitutive  constants,  u and  c.  For  the  sake  of 
brevity,  only  Le'  will  be  considered.  Once  Le'(J)  is  known,  Le"(J)  is  obtained 

AW  AA* 

by  replacing  tQ  and  uq  by  e and  u,  respectively. 

By  equation  (2-15) 

Zmn  - [ W • (JuA  + l (2-25) 

mn  j m n n 

where 

A ■ magnetic  vector  potential  due  to  J 
n n 

4 * scalar  potential  due  to  a > surface  charge 
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Applying  the  one  dimensional  divergence  theorem  to  the  vector 
that 


<t>  w , 
n m 


and  noting 


V • <(»  w - $ v • w + w 
n m Tn  m m 


the  following  relationship  is  obtained 


Define  a such  that 

SB 


• W di  - - 
m 

AAA 


W dz 

D 

AAS 


m 


(2-26) 


(2-27) 


(2-28) 


Observe  that  equation  (2-28)  has  the  fora  of  the  continuity  equation 
if  and  a are  interpreted  aa  currant  and  charge,  respectively.  An  alterna- 

u)  ID 

tive  form  for  the  [Z]  matrix  element  can  be  obtained  by  substituting  equations 
(2-27)  and  (2-28)  into  equation  (2-2S) . The  new  form  is  computationally  more 
attractive. 


Z 


mn 


Jdi 


• A + ojjdl 
n in  n 

A/wSi 


Define  the  two-dimensional  Green's  function  G(r,r'). 


(2-29) 


G(r.r')  - 

^-»<2)(k|r  - r'|) 

(2-30) 

A(r)  - u 

I J(r')  G(r,r')di* 

(2-31) 

♦(r)  - J 

| q(r')  G(r,r')dS,' 

(2-32) 

c 


whsre  q is  related  to  J by  the  equation  of  continuity 
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V ■ J - - j«q 


(2-33) 


After  aubstltutlng  equations  (2-31),  (2-32),  and  (2-33)  into  equation  (2-29), 

the  expression  for  Z becomes 
ton 

Z - j«  f { f tw  • yj  <r’)  +^qn(r,)l  G(r,r’)di'}d£ 
on  Jim  n e n 


c c 


c c 


tJwulL  • j + -i-C7.  w)(7’  • J )]G(r,r’)di,di  (2-34) 
m n i(uc  m n 

//A  A/VS  «*AA  WS  AM 


Mote  that  the  primed  symbols  refer  to  source  location  variation,  while 
the  unprimed  symbols  relate  to  variation  in  field  point  location. 

The  specific  formulation  proceeds  by  dividing  the  contour  C into  N 
segments,  not  necessarily  equal  in  length.  There  are  M segments  and  Htl  points, 


Fig.  2-1.  A cross  sectional  contour. 


r 
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and  A Is  the  path  length  proceeding  counterclockwise  around  contour  C. 

The  sets  of  expansion  and  testing  functions  are  chosen  as  triangle 
functions  for  both  electric  and  magnetic  surface  currents. 

Wk  - T(i  - ik)u4  k"l,2, . . . ,N  (2-35) 


where  Is  the  unit  vector  tangent  to  C path  length  value  1,  and  T Is  the 
triangle  function  defined  by 
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Mow,  Z can  be  written  as 
ton 


- f [JttHlW  • + T^r  (V  • W )<V*  • J )]G  dA'dl 

m j J ® “ Jwe 


Amf2  An+2 


1 [Juw  T(A-An)  T<l'-ln)(ut  • u’) 

' * AAA  AAA 


V-2  ln-2 


+ ^T'(A'  - ln)  T’(l  - *n)]Gd!*< 


(2-39) 


The  subscript  m indicates  the  nth  triangle  testing  function  and  the  eubscrlpt 
n Indicates  the  nth  triangle  expansion  function*  T'  is  the  derivative  of  the 


triangle  function. 


T' (l-»k) 


*k+2  " Ak 


Ak  " lk-2 


\ t 


1 - Ak+2 


A - Ak+2 


Ak-2  ^ A 


(2-40) 


T'(t-ik)  -i 


AAk  + AAk+l 


0 < i-lk  < 0 ik  ♦ 41^ 


AAk-2  + AAk-l 


-(At^+AA^)  < t-ik  < 0 


A“Ak  - &Ak  + AAkfl 


“<AAk-2  + Ai^)  > t-*k  (2-41) 


The  triangle  function  Is  approximated  by  four  pulses  with  amplitudes 
«nd  h^  as  shown  below 


Fig.  2-3.  A four-pulse  approximation 


where 


V-2 


y At 


k-2 


“k-2  + AVl 


fc  . 2 L\  * A*k+1 

\ Aik  + Atk+1 


‘He-! 


A*k-2  * 2 A*k-1 
Uk-2  + A*k~l 


\+l 


I A*k+1 
“k  + A*k+1 


<2- 


The  derivative  T'  of  the  Triangle  function  can  be  represented  graphically  as 


I 


AAk  + A& 


k+1 


Fig.  2-4.  Th«  derivative  of  a triangle  function. 


Consider  the  contour  Interval  spanned  by  the  expansion  or  testing 
triangle  function  as  shown  in  Fig.  2-5. 


Fig.  2-5.  Straight  line  representation  of  the  contour. 
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To  evaluate  the  Integrals  In  equation  (2-39)  each  such  portion  of  the 
contour  Is  replaced  by  straight  line  segments  drawn  between  the  points  or  the 
actual  contour  defined  by  1^,  Hj,...,  *h+1*  The  inte8r*tlon  variables  in 
equation  (2-39)  are  taken  along  these  straight  line  segments,  f or  the  nth 
expansion  function,  the  Index  p»l,2,3,4  is  associated  with  the  four  pulse 
intervals  for  increasing  path  length.  Similarly,  the  Index  q»l,2,3,4  is 
defined  for  the  mth  testing  function.  Equation  (2-39)  can  be  put  In  the 
following  form 


Z 

mn 


4 4 

l l 

q«l  p*l 


[jwW  TpTq<up 


+ -r-  T'  T']  GdJt  dl 
jut  q p p q 


(2-43) 


Note  that  Tp  and  are  the  amplitudes  of  the  pth  and  qth  pulses, 
respectively.  The  unit  vectors  up  and  u^  are  parallel  to  the  straight  lines 
of  the  pth  and  qth  Intervals  in  the  direction  of  Increasing  path  length. 

Observe  that  each  of  the  sixteen  terms  on  the  right  hand  side  of  equation 
(2-43)  results  from  one  of  two  situations.  Either  the  pth  and  qth  intervals 
coincide,  or  they  do  not.  These  two  situations  will  be  considered  separately. 


(1)  Noncoincident  Intervals 

In  this  situation  each  Integral  in  equation  (2-44)  is 
approximated  by  the  product  of  its  Integrand  evaluated  at  the 
interval  midpoint  times  the  interval  length.  Hence  equation 
(2-44)  becomes 
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4 4 


Hv2)(kR  ) 


am 


-r  -r  • »»  / 

-I  I At  At  [jun  T T (u  • u ) + T*  I']  °~  (2-44) 

«-i  P 9 P q P q jwe  p q 4j 

H ■**  P A AAA  AAA 


/SA^  AAA 


At  and  At  are  determined  by 

p q 


Atfc  - tk+1  - tfe  (k  - 1,2, 3, 4) 

and  R is  the  distance  between  the  aiidpoints  of  the  pth  and  qth 
pulses. 

(li)  Coincident  intervals 

For  coincident  pth  and  qth  intervals  the  integral  evaluations 
proceed  as  follows.  The  q integral  is  approximated  by  the  product  of 
the  Integrand,  sampled  at  the  midpoint  of  the  interval,  times  the 
interval  length.  The  p integral  is  then  evaluated  as  an  improper 
Integral. 

(2) 

The  small  argument  approximation  for  ' (kR)  is 


H<2)(kR)«^l  - j £ log 


(2-46) 


where  log  y is  Euler's  constant.  Then  for  coincident  pulse  Intervals 


Vl  Vl 


fcq  tp 


, H(2)(kS) 

tj™  y,(up  • y + jij.  i;  T-J  -Sj dy*q 


AAA  /VA 


P+1 


■ f V110"  Vi + jtr  tp  *;i  | Hi2)<kR)d* 


(2 -47) 
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Note  that  the  Integrand  la  singular  at  the  midpoint.  The  Improper  Integral 
can  be  treated  as: 


P+1 


V 2 


H<Z)(kR)d  - | tl  - J flog  <2fcM)]dx 


-At  /2 
P 


-e 


- lim 
e-K) 


[i-J:  log  (-  *&)  ]dx 


-Atp/2 
At  /2 


+ 11m  tl  - J f log  (^)]dx 

s-K)  * 1 


2 YkAt 

- 4tp[l  - 1 f log  -j^l 


(2-48) 


Therefore  Z can  be  expressed  as 
ran 


4 4 


2«  • i l I t„t>„  •»„>-  = *:  »:» 


q-l  p-1  p q p q P q 106  P q 


(2-49) 


where 


Z - 


r »'2) 


noncoincident  Intervals 


2 ykAtD 

1-J  — log  — coincident  intervals 


Equation  (2-49)  la  used  to  compute  the  two  parts  of  each  matrix  element, 
with  eq,  vjq,  and  kQ  In  one  part  and  e,w,  and  k in  the  other.  It  can  be 
readily  observed  that  the  use  of  equation  (2-49)  will  lead  to  a symmetric 
[Z]  matrix. 
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2 . 3 Evaluation  of  [El  And  [Cl  Matrix  Elements 

Matrix  elements  for  [B]  are  expressed  by  equation  (2-16) 


V - f "I  • HO  (Mn)dt 


(2-50) 


Because  of  the  discontinuity  of  the  curl  operator  at  the  boundary,  care  should 
be  exercised  in  evaluating  equation  (2-50).  The  Green's  function  is  singular, 
and  a simple  interchange  of  differentiation  and  integration  is  not  always  pos- 
sible. Note  that  the  operator  C consists  of  two  kinds  of  operators,  namely, 

C'  and  C".  C'  is  for  the  outside  field,  and  C"  for  the  inside  field,  with 
respect  to  the  material  body.  In  the  following  development,  the  symbol  C can 
be  either  C'  or  C"  unless  stated  otherwise.  Since  the  incident  field,  H*  In 
the  present  case,  is  considered  to  be  axially  directed,  there  will  be  a circum- 
ferentially directed  electric  current  and  an  axially  directed  magnetic  current. 
Let 


- id  - 


and 


(2-51) 


M ■ T(i  - £.  )u 
n k r 


Hence,  equation  (2-51)  takes  the  form 


B - 

mn 


-M 


T(£-£  )u.  • V 
tn  x 


\ 


T (£'-£  ) H<2)(kR)d£’dt 
n o 


4 4 


- -r  -T 

-■?  I l 

q«l  p-1 


T u.  • V x u f T (kR)dl'di  (2-52) 

<[  X t j p 0 


A(VS  /WA  Aae 
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Since  Tp  and  are  constant  between  tp  and  t^^,  t^  and  t respec- 
tively, hence  they  can  be  taken  outalde  the  Integral  signs,  so 


4 4 >+1  >1 

B-  ‘ ' i v,  "« 

H A P A *.  /VM  A/Vs  /w%A  a. 


4 4 

I 

q-l'  p 

*q  "p 

Again,  two  situations  will  be  treated  separately. 


(kR)di'dJt 


(2-53) 


(1)  Mon- coincident  p and  q Intervals 

The  Hankel  function  is  continuous  end  differentiable.  After 
performing  the  Indicated  curl  operation  in  equation  (2-53),  and  noting 


that  Ug  refers  to  q coordinates,  then 
4 4 3H*2)(kR  ) 


3«i2)(kR_) 


» - - T I I At  { TT-Pq  Ax  - -&—&L-  Ay  }T  T 

■“  4 q-l  p-1  p 3yq  q 3*q  q p 


(2-54) 


where 


" [(xc  - xn)2  + <yo  ‘ y„>2l1/2 

pq  q p q P 


(2-55) 


and 


3H^2)  (kR)  m x - xn 


(2-56) 


Hence 


ran 


8H‘2)(k»)  ,,,  y - y 

-*55 


.44  -IT  kH.<2)  (kR) 

i £ ,1 4S  -*V >-  % - V, 


(2-57) 


+ (x  - x )Ay  J 

p q q 


(2-58) 
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Equation  (2-58)  is  obtained  through  the  application  of  triangle 
expansion  and  testing  function  employing  a four-pulse  approximation  to 
the  triangle,  and  the  Integrand  was  evaluated  at  the  midpoint  of  each 
pulse  interval, 

(11)  Coincident  p and  q 

Note  that  in  this  case,  the  method  used  in  evaluating  the 
improper  integral  for  [Z]  can  not  be  applied  here  because  the  curl 
operator  is  not  continuous  across  the  boundary,  for  instance 


P.q ______ 

Boundary  Surface 


Fig.  2-6  Boundary  Surface 

By  visualising  a current  sheet  that  flows  into  the  paper  as  shown  in 
Fig,  2-6,  it  is  evident  that  the  tangential  field  component  will  de- 
cidedly be  sero  as  e 0. 

A better  way  is  to  find  the  field  at  a point  above  the  boundary 
surface,  then  find  the  limit  as  it  approaches  the  boundary  surface  from 
above.  Before  performing  the  limiting  process,  the  integrals  in  equation 
(2-53)  becomes 
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q+i  Vi 

| Mt  • V X uf  [ H^QtDdiMi 

' — - WV  <**A  ' 


At, 


rq+1  »H<2>  <U0  SH<2)  (kft) 

J ' 


s -]  • tu.dx,  + Ujr4y,] 


Not*  that  in  aquation  (2-58)  rafars  to  tha  q-coordinataa  and  th*  p- 
integral  is  avaluated  as  tha  product  of  tha  Intagrand  sampled  at  tha 
midpoint  of  th*  p-lntarval  tlmas  tha  interval  length. 

A local  coordinate  system  la  constructed  for  tha  evaluation  of  th* 
improper  lntagral  as  R -*•  0.  A local  coordinate  system  is  shown  In  Fig.  2-7. 


+ y2]1/2 


Fig.  2-7.  A local  coordinate  system. 
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Next,  consider  the  ;erm 


At/2  3H(2) 

J °9y 


(kR) 


dx* 


-At/ 2 


and  noting  that 


3H^2)  (kR) 


»y 


- - kH<2)(HR)  |£ 
(2) 

By  using  small  argument  approximation  for  (kR) 


H<2)(kR)  - J1(kR)  - JN1(kR) 


,M+  < 11 

2 3 r kR 


aquation  (2-59)  becomes 
At/2  3H<2)(kR) 


-At/ 2 


®y 


At/2 

dx'  - - k | 

-At/2 


rkR  . i 2 « 3R  j * 
[T  + 3 rrkR1  3y  dX 


Note  that  ei  x 0,  y -*■  0,  the  Improper  integral  approaches  -J2. 

Finally  the  expression  for  B can  be  stated  as  follows. 

mn 


i A A 

■f  I I,  VpV 

q«l  p"l 


(2-59) 


(2-60) 


(2-61) 


(2-62) 


(2-63) 


where 
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B 


B 


kH.(2)(kR  ) 

— g t-  (y_  - yJAx.  + (*_  - x Hy  ] 
Rpq  p q q p q q 


J2 


(non-coincident) 

(coincident) 


Equation  (2-63)  ie  used  to  compute  the  two  parts  of  each  [B]  matrix  element, 
R in  one  part  (due  to  primed  operator)  and  k in  the  other  (due  to  doubly 
primed  operator).  Since  C*  ie  an  outside  operator  and  C"  is  an  inside 
operator,  the  values  of  B for  ths  coincidental  case  vlll  have  opposite  signs. 
Hence  the  coincidental-pulse- Interval  situation  contributes  nothing  to  the 
values  of  the  matrix  elements. 

/V 

In  spite  of  the  fact  that  tc]  ■ [B],  it  is  advantageous  to  evaluate 

C explicitly.  The  procedures  Involved  will  be  essentially  the  same  as  in 
mn 

evaluating  B_.  Recall  equation  (2-17),  and  it  can  be  expressed  in  greater 
mn 

detail  «s 


mn 


w“  • (-jC)(J)dt 
m n 

~~  /W 


A specific  form,  suitable  for  computational  purposes  is  developed  in  a 

manner  similar  to  that  used  for  B . Considerations  governing  the  choices  of 

mn 

expansion  and  testing  functions  are  the  same  as  those  discussed  at  the  begin- 
ning of  this  section.  Note  that  the  electric  surface  current,  in  the  present 
case,  is  circumferentially  directed. 


W®  - T(* 
in 
wv* 

J - T(a 
n 

«/•*“•** 


- I.  )u 
k t 


Vu, 


(2-63) 
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where  T is  defined  by  equation  (2-38)  and  u is  the  unit  vector  in  the  direc- 

z 


tion  of  the  axis;  is  the  unit  vector  along  the  cross-sectional  contour.  C( 

can  be  expressed  as 


mn 


on 


, 4 4 V1 

H q-1  p-1  J 


T u • V « 
q z 


“J 


P+1 


T H*2j(kR)dJl’d8, 
P 0 


(2-66) 


q P 

The  evaluation  of  the  Integrals  appearing  in  equation  (2-66)  is  facilitated 

by  approximating  the  triangle  function  by  four  pulses.  The  index  p ■ 1,2, 3,4  is 

associated  with  each  pulse,  respectively,  for  the  nth  expansion  function,  while 

the  index  q - 1,2, 3,4  is  similarly  defined  for  the  mth  testing  function.  Since 

T is  constant  between  t and  t . , and  T is  constant  between  t and  t , they 
P P P+1  q q q+l 

can  be  taken  outside  the  Integral  signs.  is  the  unit  vector  along  the  contour 
with  respect  to  p coordinates.  Hence 

t 


.44  (^1 

C»n  ' - I l,  l I W 

q-1  p—l  ’ • ^ r 


(2) 


(kR) 


,44  r 3H’  • (KRJ 

- - T l l T T At  t— dy  - 

4 q-1  p-1  p q * > 3xq  p 


+ udypH<2)(kR)) 

'P*1  3H*2\kR)  3H*2)(kR) 

” o 


3y, 


dx J (2-67) 

P 


Hots  that  the  integral 


P*1  3H(2)(kR)  3H^2^  (kR) 

l — It: dy  - — dx  ] 

3x  p 3y  p 


(2-68) 


can  be  evaluated  just  liks  in  the  previous  B esse.  It  follows  that  equation 

on 

(2-68)  becomes 
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-“i”  ^ 


- X. 


-9-  Ay  + 


pq 


y„  - y. 


pq 


v 


(2-69) 


For  coincident  p~pulse  and  q-pulaa  intervals,  the  evaluation  of  the  Improper 

lntegrel  is  Identical  to  that  developed  for  B . To  this  point,  the  matrix 

nm 

elements  for  [C]  can  be  conveniently  specified  as 


C 

an 


* l 


q.l  p 


A 

I At 
1 


T T C 

q p q 


(2-70) 


where 


tu(2)  IkRl 

c \ ~ [-  (x  - x )Ay  + (y  - y )Ax  ] 

p q p p q p 


pq 


(noncoincident  Intervale) 


C - -j2  (coincident  intervela) 


Equation  (2-70)  is  used  to  compute  the  two  parts  of  each  [C]  matrix 
elements;  one  part  is  due  to  the  outside  operator  C' , and  the  other  Is  due 
to  the  inside  operator  C".  Again,  for  the  coincidental-puls e-interval  situ- 
ation, the  net  contribution,  to  the  value  of  each  matrix  element,  is  sero, 
because  the  two  values  of  C in  equation  (2-70)  have  opposite  signs, 

2.4  Evaluation  of  [Yl  Matrix  Elements 

The  only  expression  left  to  be  developed  1b  that  for  the  [Y]  matrix 
elements.  By  equation  (2-18) 

Y -[  VJ°  • L (M  )di 
nm  I m an 

- f w®  • (jwF  + V4®)dA 

] m n n 

AM  MA  AM 


(2-71) 


The  superscript,  m,  Indicates  magnetic  quantities.  By  equations  (2-26) 
to  (2-28),  the  following  ia  obtained 

Y “ Ju  f (.V  • F + oJ>*)dl  (2-; 

nm  j m a m < n 

In  this  case  all  the  currents  are  axially  directed.  The  expansion  and 
testing  functions  are  chosen  as 


wk  - ^ - T(i  - ak)ug  (2-: 

AM  /AAS  AAA 

where  T is  the  triangle  function  defined  by  equation  (2-38),  and  u is  the 

s 

AM 

unit  vector  in  the  axial  direction.  The  continuity  equation  in  this  case 
is 


a 


n 


M 

n 

AM 


(2-: 


Note  that  Mn  is  u,  directed,  so 
XI  s 

AM  AM 

v . Mn  - 0 (2-; 

and  it  follows  that 


■ 0 and  $m  ■ 0 (2-1 

n n 

Therefore 

Y - Jw  f w“  • A®  di 

mn  inn 

q <**■*  AAA 

“T  f f <wn  • Mn)H^2)(kR)di’di  (2-; 

' ' (W,  A.--S 

c c 


where  the  unpriaed  integration  la  taken  over  field  points  and  the  primed 
integration  over  the  source  points.  A specific  form  is  developed  in  a 
manner  similar  to  that  used  for  Z . Equation  (2-77)  can  be  exprassed 
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In  greater  detail  as 

*mf2  *n+2 

Ytan  " T j { T(A'-r)  H^OOOJdi’dA  (2-78) 

*m-2  *n-2 


The  evaluation  of  the  integrals  appearing  in  equation  (2-78)  ia  carriad 
out  by  approximating  the  triangle  function  by  four  pulses  as  ahown  in  Fig.  2-3. 
Equation  (2-78)  can  be  written  as 


4 4 >l  V1 

’.-f  uj 

q-1  p**l  i > 


T T H^OclOdA'dl 
p c o 


(2-79) 


The  indices  p and  q have  the  usual  meaning.  T and  T have  already  been 

P <1 

defined  by  equation  (2-42).  Each  of  the  sixteen  terms  contributing  to  Y 

mn 

fells  into  one  of  two  categories.  Either  the  p-pulse  end  q-pulse  intervals 
coincide,  or  they  do  not.  If  the  latter  is  true,  each  integral  is  approxi- 
mated by  the  product  of  its  Integrand  sampled  at  the  interval  midpoint  times 
the  interval  length.  The  expression  for  Y becomes 


mn 


- *r  t J T T At  At  HW(kR  „) 

4 q-1  pSl  P « P *»  0 W 


(2-80) 


As  previously  stated,  Atp  and  At^  are  defined  by  equation  (2-45)  end 
Rpq  la  the  distance  between  the  midpoints  of  the  p-pulsa  and  q-pulse  intervals. 
For  coincident  p-pulsa  and  q-pulse  intervals,  the  improper  integral  and  its 


evaluation  era  the  same  as  those  in  Z 


mn 


Hence 
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Vi  Vi 


TT  H^(kR)di'di 
p q o 


c t 

q p 


1 Yk  4 t 

Vq  V\  U " J f 108  ("^r£)1 


(2-81) 


where  logy  - Euler's  constsnt 
e * 2.7 18  .... 


The  flnel  expression  for  Y Is 

on 


where 


4 4 

r - x y y tt  ititu 

m 4 q-l  p-1  P « P « 


Y ■ H^(kR  ) (non-coincident  intervals) 
o pq 

, Y^4t 

» [1  - J ^ log  (■y-f)  ] (coincident  intervals) 


(2-82) 


2.5  Excitation  Matrix.  Measurement  Matrix,  and  Scattering  Cross  Section 

The  matrix  elements  of  the  excitation  matrix  ere  represented  by  two 
expressions,  equations  (2-19)  and  (2-20).  It  is  Important  to  reallsa  that 
the  transformation  of  (2-19)  and  (2-20)  into  computable  forms  depends  on 
the  type  and  polarization  of  the  impressed  field.  In  the  case  under  con- 
sideration, the  excitation  is  assumed  to  be  a z-dlrected  magnetic  field 
of  unit  magnitude.  The  incident  field  is  given  by 


H^p)  - e’Jk**&  (2-83) 

z 

The  wave  number  vector  k points  in  the  direction  of  travel  of  the  Incident 
wave.  A coordinate  system  for  the  evaluation  of  the  excitation  matrix 
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elements  Is  shown  In  Fig.  2-8. 


Fig.  2-8.  Incident  field. 


Equation  (2-19)  will  be  considered  first.  The  testing  function  is 

«£  - T<a  - *k>UA  (2-84) 

For  plena  wave  excitation  the  ^-directed  electric  field,  associated  with  the 
z-dlrected  magnetic  field  defined  by  equation  (2-83),  is 


K*  - ~ riu.e 

11 

± i ^k(xmBC0,  ♦ + ymoBln  * > 

■ - n(-  u sin  $ + u cos  Oe  **  ** 


(2-85) 


where 


u.  * ^-directed  unit  vector 
ux  ■ x-directed  unit  vector 


y-directed  unit  vector 


i! 


$ * incident  angle 

x ,y  ■ midpoint  coordinates  of  each  straight  line  segment 
mp  mp 

ri  - intrinsic  wave  impedance 

By  using  four-pulse  approximations  for  the  triangle  testing  functions, 
equation  (2-19)  can  be  expressed  as 


vi-  I 


ua  T • E (p)d& 

AAA  Y 


(2-86) 


Mote  that  the  excitation  matrix  element  V is  given  as  the  component  of 

m 

E1(p)  tangent  to  the  contour  for  the  nth  triangle.  The  integral  in  equa- 
tion (2-86)  is  evaluated  as  the  tangential  field  component  of  E*  sampled 
at  the  midpoint  of  each  p-pulse  interval.  Hence 


vi  - - n l T { 

“ p-1  p 

jk(x  cos  + y sin  $*)  . . 

e ”p  mp  t-  Ax  sin  + Ay  cos  *4]  (2-87) 


where  Ax  and  Ay  are  the  rectangular  components  of  the  pulse  Interval.  A 
portion  of  the  contour  is  shown  in  Fig.  2-9  which  illustrates  how  equation 
(2-87)  is  obtained. 

Equation  (2-20)  can  be  evaluated  in  a similar  manner.  The  testing 
functions  ere  triangle  functions,  and  each  triangle  function  is  represented 
by  equation  (2-39)  and  chosen  to  be  z-directed. 


\ - T(£  - *k)uz 

.vAA  AAA 


(2-88) 
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mp 


Fig.  2-9.  A partial  contour. 


The  evaluation  of  (2-20)  Is  quite  straightforward  and  the  procedures  are 
identical  to  those  used  In  the  evaluation  of  V*.  Hence 

IQ 


. 4 jk(::  cos  + y sin  (j)1) 

I1  - j 7 T At  - “P  “P 

m L - « r 


P-1 


p p 


(2-89) 


The  distant  scattered  field  can  be  evaluated  by  reciprocity.  A z-directed 


magnetic  current  filament  at  p of  strength  M is  adjusted  to  produce  the  unit 

o P 

plane  wave  Incident  on  the  material  body 


. jk(x  cos  $8  + y sin  <t>8) 

H1  - u e n n 

Note  that  M produces  a ^-directed  E1  and  a z-directed  H* 


(2-90) 


:\ 


Fig.  2-10.  A two-dimensional  contour  and  z-dlrected 
magnetic  current  filament. 


By  reciprocity 


or 


where 


Equation  (2-7) 


it  la  evident  that 


“ Hn  * if-  [ (E1  • J - H1  • M)dl 
p n J e/A  /va  /*A 


-«c-r 


r\f 

f 


dl 


£ - 


J 


jM 


can  be  expressed  in  matrix  form  as 


f - 


IV J tMj 
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With  Che  help  of  equation  (2-93) , a new  form  for  equation  (2-92)  la 


p 


(2-94) 


Note  that  [I  ] and  [V  ] are  column  matrices,  and  the  matrix  [D]  la 
n n 


ID]  - [[(  eJ  • [Jjdl]  [f  [M^l  • jH^dtn 

r ***  r 


(2-95) 


The  values  of  I and  V can  be  obtained  from  equation  (2-14)  by  matrix 
n n 

Inversion.  The  constant  1/M  Is  that  needed  to  produce  a plane  wave  of  unit 

P 

amplitude  at  the  origin,  which  la 


_1_ 

M 


■fern?*  too’ 
-f 


(2-96) 


Redefine  equation  (2-95)  as 


where 


[D]  - ([De] 

[d“]> 

_e 
D - 

( E1  • 

[Jjdt 

n 

C 

n 

_m 

D ■ 

(jH* 

• IMJdl 

n 

J /M 
c 

n 

ASS' 

(2-97) 


(2-98) 


(2-99) 


The  evaluation  of  the  Integrals  In  equation  (2-98)  and  (2-99)  Is 
straightforward,  and  the  procedures  Involved  are  completely  analogous  to 
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those  used  in  the  evaluation  of  V*  end  I*.  Only  the  results  will  be  given 

m m 


here 


e ft  Jk(x  cos  * + y««8in  ♦'  > 

D*  - - h I T e np  np 

n p-1  P 


_ $ Jh(x__co*  + y__«in  <>*) 

D'-j  Tit.  °P  np 

n p-1  P P 

The  scattered  field  can  be  axpr-rred  as 


[-Ax  ein  $ + Ay  cos  A ] (2-100) 

(2-101) 


-f  h<2)0»o)[[dJ]I^)1 

‘[I  J 

‘ n' 

or 

[V  1 
1 nJ 

• m 

"[Z] 

IB]" 

-1 

’<<> 

[C] 

m 

[V] 

(2-102) 


(2-103) 


In  the  scattering  problem,  the  blatatic  scattering  crose  section  a is 
a parameter  of  interest.  It  Is  defined  as  the  width  for  which  the  Incident 
wave  carries  sufficient  power  to  produce  the  field  Ep,  Hp  by  omnidirectional 
radiation.  It  may  be  expressed  as 


or 


o(<t>8)  - 2ir  p0 
lim  p„ 


■ 2 tr  p0 
Urn  po+  • 


Ep(»") 


Hp(^) 


(2-104) 


;\ 
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i 

i 


(2} 

The  large  argument  approximation  for  H (kR)  la 
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CHAPTER  3 

CHARACTERISTIC  MODES  - A SURFACE  FORMULATION 
3.1  Theoretical  Development 

The  treatment  of  characteristic  nodes  for  perfectly  conducting  bodies, 
starting  from  the  Impedance  operator  for  the  conducting  surface,  has  been  given 
by  Harrington  and  Mauts  [2].  In  terms  of  the  polaritetlon  current  and  the 
magnetization  current,  a volume  formulation  of  the  characteristic  mode  theory 
for  dielectric  and  magnetic  bodies  has  also  been  treated  [4].  In  this  chapter 
a theory  of  characteristic  modes  for  material  bodies  (dielectric,  magnetic,  or 
both)  based  on  a surface  formulation  is  developed.  The  appropriate  operator 
formulation  of  the  problem  is 


(3-1) 


To  emphasize  the  eynmetric  nature  of  the  matrix  of  operators,  the  off-diagonal 
operators  are  denoted  by  a single  symbol,  N.  Define  the  following  rlee  vectors 


- - 

- - 

j 

Le  N 

j 

E1 

MA 

m 

/V* 

1 

a 

55 

A 

JH1 

LV*  j 

f - 


and  the  matrix  of  operators 


& 


S ■ 

<*/*> 


E 

aaA 


JH 


(3-2) 


T - 


Le  N 


N Lm 


- -JC. 


(3-3) 


4 


where  N 
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Equation  (3-1)  can  then  be  written  as 


Tf  - g 
^ (*//■> 


Define  the  symmetric  product 


<f.  & ■ 


fg  da 


■fl 


(J  • E - M • H)ds 

A***  ^ 


C3-4) 


which,  for  f a source  quantity  and  g a field  quantity,  la  reaction.  The 
product 


<£* 


,»>  - ( [ 7*g  da 

A*  J J A*. 


'll 


(J*  • E + M*  • H)ds 

AM  VA  /yw\  vV’ 


(3-5) 


ia  a aultable  inner  product  for  the  Hilbert  apace  of  functions  f,  g in  S. 

aaA 

If  f la  a source  quantity  and  g a field  quantity,  the  real  part  of  (3-5)  is 

A"* 

time  average  power,  but  the  Imaginary  part  of  (3-5)  differs  from  the  usual 
Imaginary  power.  It  ia  easy  to  show  that  T ia  symmetric,  that  la, 


<f^,  Tf 2>  ■ <fj,  Tf^>  by  reciprocity.  The  operator  T can  be  expressed  in 

AAA  AAA  AV  AAA« 

terms  of  its  Hermit ian  parts  as  T « T^  + JTg  where 


Tx  - j (T  + T*)  - 


t2  ■ Ji  <T  ■ T*} 


R 


N1 

G 


(3-6) 


(3-7) 
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Here  and  Nj  are  the  Hermitian  parts  of  N,  R and  X are  the  Hcrmltian  parts 
of  Z,  G and  B are  the  Hermitian  parte  of  Y. 

By  aquation  (1-9),  the  fields  due  to  J and  M can  be  expressed  as 


- 

m - 

— 

J 

|e  1 

Le'  -JC* 

L’ 

.1M 

■ mm 

5 l1 

-JC  La' 

(3-8) 


Aa  far  as  radiation  is  concerned,  the  contribution  due  to  the  doubly  primed 
operators  is  sero.  The  power  radiated  by  any  J and  M on  S is  given  by 


K»<P,> 


(E  • J*  + H * H*)d« 

AAA  A** 

(E  * J*  + M*  • H)ds 

MM  AAA  AAA  AAA 


- Re  (<f*,Tf>)  (3-9) 

AAA  AAA 

Hence  the  time  average  power  delivered  by  a aource  f is 


Re(P  ) - Re  <f*,Tf>  (3-10) 

• AAA  aAA 

The  imaginary  part  of  <f*,Tf>  is  not  simply  related  to  reactive  power. 

AAA  a AA 


Using  six-vector  notation,  we  formulate  a theory  of  characteristic 
modes  which  parallels  that  of  the  volume  formulation  [4].  The  eigenvalue 
equation  defining  the  modes  is 


T,(f  ) - X T.(f  ) (3-11) 

& n n x n 

where  and  T2  are  real  symmetric  operators.  Hence,  all  eigenvalues 

are  real  and  all  characteristic  sources  f may  be  chosen  real.  In  expanded 

n 


form 
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f ■ 
n 


J 

n 


JMfl 

L»~\J 


(3-12) 


which,  for  character istlc  sources,  implies  that  N is  imaginary  and  J is 

n n 

AAA  /VVN 

real.  The  characteristic  sources  can  he  normalized  to  radiate  unit  power, 
and  the  usual  orthogonality  relationships  expressed  as 


<f*. 

T,f  > 

- <f  , T.f  > 

m 

1 n 

m*  In 

*** 

***> 

<f*. 

T,f  > 

- <f  , T„f  > 

m 

2 n 

nr  2 n 

AM 

/VS^  /v* 

<f*. 

Tf  > ■ 

■ <f  , Tf  > ■ 

m 

n 

m n 

f 

ran 


n mn 


n mn 


(3-13) 


where  6 is  the  Kronecker  delta.  The  field 

Ou 


IK 


(3-14) 


due  to  a source  f is  called  a characteristic  field.  In  the  radiation  zone 
n 

AM 

the  charactsrlstic  field  is  of  the  form  of  an  outward  traveling  wave,  and  it 
is  completely  characterized  by  either  E or  H . 

^ A«J» 

Let  f and  f be  two  characteristic  sources.  By  equation  (3-13),  the 

u n 

AvSA 


Le  N 
N Lm  jM_ 


following  expression  la  true. 


<f  , Tf  > -<  [J  JM  ] 
on  n n 


* 0 for  ra  ^ n 

Equation  (V13)  is  essentially 

f (J  • E - H • M )ds  ■ 0 
JJ  m n n m 


(3-15) 


0-36) 


/\M 
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where  E and  H are  produced  by  f . Because  J is  real  and  M is  imaginary, 
n n n m nj 

wa  have 

7 <J*.En  + H .M*)ds-0  (3-17) 

j j ran  n m 

A«vs  /W*  A/V*  /W* 

It  follows  that 

Ra  |f  (En‘Jm  + Hn’Vd*  * 0 (3-18) 

*vV\  /vv\  /W->  AA*i 

which  naans  that  the  reel  part  of  the  cross  power  is  zero.  In  tha 
radiation  zone  the  characteristic  waves  are  of  tha  font  of  outward 
traveling  wave,  l.e. 


E - n H x n (3-19) 

n n 

AAA  e^A  AAA 

where  n is  the  unit  radial  vector  on  SM,  The  real  part  of  the  cross  power 

AAA 

can  be  expressed  as 


80  ||  Em  * Hn  ‘ d*  " R*  | Imfi.  m ® 

/vw  ^ 


(3-20) 


The  real  part  of  tha 
zero.  Hence, 


cross  power  between  j(E  ,H  ) and  (E _,H  ) is  also 

n n n n 

AAA  AAA 


aAa 


AAA 


ds  * ri 


AAA 


“nd* 


6 


nm 


(3-21) 


3.2  Characteristic  Equation  and  Modal  Representation 

In  the  preceding  section,  the  analytical  development  was  based  on 
the  interpretation  of  operators.  The  reduction  of  operator  equations  to 
matrix  equations  can  be  effected  in  tha  usual  manner  by  the  method  of 


moments. 
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Ut 


fn  * 

AAA  J 


< y;  ♦ vjtj  ) 


where 


n- 


wj! 


aAA 


-m 

fJ  " 

VA 


After  substituting  equation  (3-22)  into  equation  (3-11) , the  following 
ia  obtained. 


Perform  inner  product  with  electric  testing  function  VfJ, 

AAA 

« i v*  vj  * + 1 v "?•  vj  > > 

j AAA  j AAA.  AAA 

- *„  < I y «J.  vj  > + J < wj.y*  > ) 

J A/vs  AA  . J AAA  AAA 

and  with  magnetic  testing  function  V^. 

ASA 

( Jlj  < »J,  IjtJ  > + Jvj  < Vlj,  yj  . ) 

"*  aAA  AAA  ^ AAA  ASA 

- x„  ( £ y b“.  vj>  + I VJ  «>?.  Tlfj  » > 

J AAA  AAA  J AAA  MVS 


Equation  (3-25)  and  equation  (3-26)  can  be  put  into  one  matrix  equation. 


'[X] 

[■*1 

til" 

- x 

'[R] 

'll) 

n 

K> 

181 . 

m 

n 

- 

IG] 

m 

(3-22) 


(3-23) 


(3-24) 


(3-25) 


(3-26) 


(3-27) 


< 


\ 
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The  definitions  of  [X],  [N2],  [B],  lR],  [Nj,  [G],  [I],  and  [V]  are 
obvious  by  comparing  equation  (3-27)  with  aquations  (3-25)  and  (3-26). 

Equation  (3-27)  is  the  eigenvalue  aquation  which  will  be  used 
in  the  actual  computation  of  the  modes.  In  abbreviated  form,  it  becomes 

I I2  1 t f„  1 - Xn  t TX  ] [ ] 

AAA  AAA 

Now | with  the  understanding  that  X and  fn  can  be  found 9 the  modal 
solution  for  f can  be  expressed  as 

AAA 

£ " l 

/vv»  n 

Recall  that 


Tf  - s 

After  substituting  equation  (3-29)  into  equation  (3-30)  and  performing 
the  inner  product  with  fm,  the  following  equation  results. 

I “»  * '«•  Tfn  * ' < V ’ 

Apply  the  orthogonality  relationships  given  in  equation  (3-13) . It 


follows  that 


< *..  8 > 


(3.  + J X ) < f , Tf  > 
n — » 


Explicitly, 


< fm-  " ■ I h < fI-  ” + I Vi  * *?•  » 

/vA*»  ^AA  • Aaa  /AA  X > eA  AAA 


AiV-  ^ 


The  matrix  equivalents  of  the  orthogonality  relationships 
for  the  characteristic  currents,  equation  (3-13),  ere  also  of  interest. 


Tor  example,  that  for  T^  is 


K 
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‘ V V«  * ■ * l ' h'l  + >»•  T1  ? < V*  + V?  >n  > 

W /w*  -1  M'  ^ 

* H * Vj  * *i>  Vj  - + Vi  ' £I-  Vj  > 

1 J ***  ^v\  /w<.  /w\ 

+ Vj  * fr  Vj  > - Vj  - *i«J  > > 


W AAA 


AAA  AA* 


I I ]„t  * 1 I x ],+  I I )„[  »!  1 [ V 1 

+ i 1 1„+  (Vy  o 1 1 v i 


n 


I fm  It  T1  II  £n  1 - 5mn 

AAa 


where  denotes  transpose.  Similar  derivations  hold  for  Tg  and  T. 


(3-34) 


3.3  Linear  Measurement 

Any  scalar  p linearly  related  to  the  generalized  current,  l.s. 
a linear  functional  of  the  equivalent  electric  and  magnetic  currants,  will 
be  called  a linear  measurement  of  the  current. 

Any  linear  functional  of  f can  be  expressed  as 

P ■ < g“,  f > (3-35) 

Aas  /VA 

where  gm  is  a vector  function  which  consists  of  an  electric  field  and  a 
magnetic  field.  By  equations  (3-32)  and  (3-33),  the  linear  measurement 
of  f can  be  stated  as 

< f , g*  > 

- I caftljmc < g®,  f > 

n (1  + jXn)<  fn,  Txfn  > Ji. 

AAA  AAA 


P 


(3-36) 
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where 


< g“.  ffl  > ■ I I±  < g®,  £®  > “ I vi  < s“*  > 


A«>  11 

and  define  the  following 


modal  measurement  coefficient 


„m  o , 

Kn  ■ <«  • V 

iW\ 

<gi,  f > * modal  excitation  coefficient 

n 

Equation  (3-36)  ia  a symmetric  bilinear  functional  of  g1  (the 
impressed  field)  and  of  gm  (the  measured  field).  The  symmetry  of  (3-36) 
is  a confluence  of  the  symmetry  of  the  original  operator  T.  Equation 
(3-36)  can  be  expressed  as 


I 

n 


K1  Km 
n n 


1 + J A. 


,.i 


similarly,  in  terms  of  equation  (3-32)  becomes 


(3-37) 

(3-38) 

(3-39) 


(3-40) 


1 + j \ 


n 


and  equation  (3-29)  will  take  the  form 


f - l 7 - ■ - f 

_ *•  1 + J X n 
/w-  n n ^ 


(3-41) 


(3-42) 


3.4  Characteristic  Fields  and  Scattering  Cross  Section 

The  characteristic  fields  are  linearly  related  to  the 
characteristic  currents,  fn,  and  hence  can  also  be  expressed  in  modal  form. 


B " £ 1 + j X *n 

ft  ft  /vsA* 


(3-43) 
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When  K and  f are  known,  the  field  pattern  can  be  obtained  by 
n n 

employing  equation  (3-43).  A convenient  way  la  to  evaluate  the  modal 
measurement  coefficient  first.  In  the  two-dimensional  case,  consider  a 
magnetic  current  filament,  M ■ M uffl  «t  ( p,$  )on  S^.  See  Fig.  3-1  below. 


Fig.  3-1.  A coordinate  system  for  modal  measurement 


coefficient 


By  reciprocity,  it  is  readily  seen  that 

- H . M - ( J . E - M . H ) ds  (3- 

n n m n m 

^ ^ AA  Wi 

wiiare  Hn  is  the  characteristic  field  components  produced  by  the  mode 

current  f . E and  H are  the  fields  due  to  the  magnetic  current  filament, 
n in  in 

AAA  AM  MA 

M.  To  simplify  the  analysis,  the  magnitude  of  the  magnetic  current,  M,  is 

MA  /VA 

adjusted  to  produce  a plane  wave  on  the  material  body,  i.e. 

- j k . r 

H ■ u e (3* 

tn  tn 

N*\  AAA 


Em  “ n Hm  * “km 


(3-46) 
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where  n is  the  wave  impedance  and  is  the  unit  vector  in  Che  direction  of 


/N~S 


propagation.  The  right  hand  aide  of  equation  (3-44),  in  matrix  form,  is  the 
modal  measurement  coefficient.  Hence, 

K®  - f ( J .E  -M  .H  ) 
nj'n  m n m 


il 


/AAA  AAA  AAA 


(3-47) 


Explicitly,  the  electric  field  and  the  magnetic  field  can  be 
extracted  from  equation  (3-43)  as 


E - l 

aaa  n 


n 


1 + j A n 
n /wa 


(3-48) 


H - l 
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H_ 


1 + j “n 
n 


(3-49) 


Since  the  magnetic  field  is  currently  under  consideration,  only  equation 
(3-49)  will  be  used.  The  component  of  the  magnetic  field  on  u is 

Dl 


H • um  " l 
aaa  a*a  n 
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i + J \ 


i K« 

1 r n n 

M S 1 + ^n 


. -M. 


H 


(2) 


Ocp)  l 
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u 


K1  K" 
_ n n 

1 + J * 


n 


(3-50) 


Note  that  K*  is  of  the  same  functional  form  as  K® . Equation  (2-98) 
n n 

m 

has  been  used  in  deriving  equation  (3-50). 

A commonly  used  parameter  in  plane  wave  scattering  problems 
is  the  echo  arse.  In  two-dimensional  problems  the  quantity  " echo  width  " 
corresponds  to  the  " echo  area  " of  the  three-dimensional  problems. 
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The  echo  width  is  defined  in  equation  (2-106). 

a “ 2 ir  p | H . u | 

“ w*  ® 
limit  p • 


(3-51) 


By  equations  (3-50)  and  (2-106),  the  following  expression  for  the  scattering 
cross  section  Is  obtained. 


K1  K° 
n n 


4 n4 


I T S_ 

4 1 + i 


(3-52) 
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3.5  Computational  Considerations 

The  solution  of  the  matrix  eigenvalue  problem,  equation  (3-28), 
will  be  discussed. 

[ T,  ] [ £ 1 - X [ T.  ] t f ] (3-53) 

Note  that  the  subscript  n has  been  dropped  for  brevity.  The  conventional 
method  for  reducing  (3-53)  to  a symmetric  unweighted  eigenvalue  equation 
requires  [T^]  to  be  positive  definite.  In  theory  [T^]  is  positive  aemi- 
deflnlte,  but  because  of  numerical  inaccuracies  it  is  actually  indefinite, 
with  some  small  negative  eigenvalues.  If  the  values  of  the  matrix  elements 
cover  a very  wide  range,  scaling  will  become  desirable.  The  magnitude  of 
the  scale  factor  can  be  chosen  as  such  that  all  scaled  matrix  elements 
will  be  brought,  as  close  as  possible,  to  the  6ame  order  of  magnitude. 

The  conventional  method  will  be  modified  as  follows. 

Let  [D]  be  a diagonal  matrix.  After  premultiplying  by  [D], 


equation  (3-47)  becomes 
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[ D ] [ T,  ] [ f ] - \[  D ] [ T.  ] [ f ] 
/w^  1 />**• 


Then  observe  that 

[ D ]t  T,  ][  D ] ( [ D 3-1[  f ] ) 

4 /W' 

- X [ D ][  T,  ]l  D ] ( t D ]-1[  f ] ) C 

The  eigenvalue  equation  aa  given  by  (3-55)  will  have  the  aeme  eigenvalues 
ae  the  original  unsealed  equation,  but  the  eigenvectors  will  be  different# 
In  other  words  the  eigenvalues  are  not  affected  by  the  diagonal 
transformation.  The  original  eigenvectors  will  be  modified  by  [D]  inverse. 
If  the  scale  factor  is  s,  [D]  can  be  chosen  as 


t D ] 


1/ 8 0 

as 

o 1 


By  equations  (3-27),  (3-55),  and  (3-56),  the  scaled  eigenvalue 
equation  is 

t T?  ][  f*  ] - A[  T*  ][  f8  1 


[ Tj  ] 


( tJ  ] - 


[ x ]/s2 

t N2]/sj 

t^  1/a 

t B ] J 

t R ]/b2 

[ ]/. 

3/s 

t G ] 

(3-59) 
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t f8  ] - 


at  I 1 
t V ] 


(3-60) 


Note  that  [I]  and  [V]  ( J and  1M)  should  first  be  recovered  from  the 

MVS 

scaled  eigenvectors  before  computing  the  surface  currents  and  the  scattered 
fields. 

Rewrite  equation  (3-57)  below 


[ T,  ][  f ] - Xt  T.  H f ] (3-61) 

* ***  1 AAO 

Note  that  the  superscripts  have  been  dropped  for  brevity.  An  approximation 
will  be  used  in  finding  the  eigenvalues  and  eigenvectors.  The  eigenvalue 
equation 


[ I,  3 r ■ y r 

* AM 


(3-62) 


is  used  to  find  a set  of  basis  functions  fot  the  vector  space.  An 

orthonormal  set  of  vectors  can  be  obtained  by  using  the  vectors  { r^  }. 

Let  { } be  the  set  of  orthonormal  vectors,  and  let  [U]  be  the  orthogonal 


matrix  which  diagonalizes  [T^]  according  to 


tTj[  Tx  ][  U 3 


U,0  0 0 

oVo  0 
o oVo 

0 0 0 u4 


(3-63) 


where  the  are  the  eigenvalues  of  [T^J  ordered  >,  >.  Wj  >,  >.• . 

Every  column  of  [U]  is  in  { U^} . Only  the  larger  can  be  considered 

accurate.  All  y^  > My^  are  put  in  [ y^  ] where  M is  some  small  positive 

number  set  by  the  estimated  accuracy  of  [T^],  Usually  M is  anywhere 
—3  —6 

between  10  and  10  . The  diagonal  matrix  (w3  i«  than  partitioned  as 
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where 


[ v 


I ux  1 [ 0 ] 

[ o ] [ u2  1 


( U5 


nrf2 


Now  consider 


n 

f - I x 
A-l  x 


U 


(3-64) 


(3-65) 


(3-66) 


(3-67) 


where  Is  e column  vector  of  [U].  This  is  a valid  expansion  because 

the  { Uj,  } vectors  form  a basis  for  vector  space.  In  matrix  form  equation 

(AN 

(3-67)  becomes 


f - [ U ] [ x ] (3-68) 

/SNA  AAA 

If  [ Uj  1 1b  set  to  zero,  it  follows  that  all  column  vectors  of  [U] 
corresponding  to  all  e [ Uj  1 *re  ln  the  nuH  “P*ce  o£  Thi8  is 
illustrated  in  Fig.  3-2. 

The  expression  for  f as  given  in  equation  (3-67)  can  be  written 


i 

i 


i 


as 


m n 

l + l Vk 

i-1  k-nrt-1 


(3-69) 


4 
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Fig.  3-2.  vector  space 


The  column  vector  [x]  in  equation  (3-68)  can  be  partitioned  as 
^ r 

t xx  ] 


t x ] - [ U ][  f ] 


[ x2  ] 


(3-70) 


where  [x^]  and  [x2l  are  column  vectors,  and  they  are  obtained  by 

<V»' 

partitioning  [x]  according  to  equation  (3-69).  Pramultlply  equation  (3-61) 


by  [U]  and  use  equations  (3-63)  and  (3-68).  The  result  is 
[Tit  T,  ][  U )[  x ] - Xt  p ][  x ] 

AAA  <S»A 

Set  [ p2  3 equal  to  zero  and  partition  all  other  matrices 
conformably.  The  following  two  matrix  equations  are  obtained. 

[ It  Xj  | + i Ay  II  *2  1 ■ M lij  ][  x1  ] 


I Ai2  It  Xj  1 + [ Ajj  1[  ^ ] > 0 

aAA. 


Note  that  [ A ] ■ [ U ][  T2  ](  U ].  Equation  (3-73)  can  be  solved  for 
[ x2  ] and  the  result  substituted  into  aquation  (3-72)  to  get 

***** 

l “ ^12^22^12  ^ 3 * At  Pi  ][x1  ] 


‘i 


'1  “1 
AAV- 


The  brackets  of  submatrices  have  been  dropped  to  conserve  space. 


(3-71) 


(3-72) 

(3-73) 


(3-74) 


•A' 

,;is 

4 
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Now  [v^]  has  only  positive  diagonal  elements  as  defined  by  equation  (3—65) . 
Observe  that 


[ ] * l yj/2]  t wj/2] 


(3-75) 


where 


l u^l 


1/2 


1/2 


1/2 

w3 


M1/2 

pm 


(3-76) 


By  equations  (3-74)  and  (3-75)  a new  and  unweighted  eigenvalue 
equation  la  obtained. 


where 


I B ][  y ) - X[  y ) 

(3-77) 

[ y ] - [ u}/2][  x,  ] 

(3-78) 

[ B 3 ■ I - A^2^22A12  ^ ^ U1  J 

(3-79) 

The  eigenvalues  of  equation  (3-77)  are  the  smaller  eigenvalues  of 
the  original  equation  (3-61),  and  the  eigenvectors  of  (3-77)  give  the 
corresponding  eigenvectors  of  equation  (3-61)  according  to 


f - t U ](  x ] - [ U ) 


( * 3 

t"A22A123 


I w:1/2][  y ] 


(3-80) 


where  [ 6 ] Is  the  identity  matrix. 


1 


Once  the  eigenvalues  and  the  eigencurrenta  are  known,  the 
equivalent  surface  currents  and  scattered  fields  can  be  obtained  by 
employing  appropriate  formulas  for  those  quantities. 
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CHAPTER  4 
RESULTS 

The  results  of  fsr  field  scattering  calculations  for  some  material 
cylinders  are  presented  In  this  chapter.  Equations  used  are  those  developed 
In  Chapter  2 and  Chapter  3. 

The  far  field  scattering  patterns  of  circular  material  cylinders  have 
been  computed  end  the  results  are  shown  In  Figures  4-1  through  4-16,  for 
perpendicular  polarization  (TE) . Figures  4-17  through  4-22  give  the  results 
for  parallel  polarization  (TM).  All  results  are  compared  with  exact  harmonic 
series  solutions  [7].  Figures  4-1  to  4-5  are  obtained  by  using  15  triangle 
expansion  functions.  Twenty  expansion  functions  have  been  used  in  obtaining 
Figures  4-6  to  4-22.  In  all  figures  the  computed  scattering  cross  section 
are  normalized  by  na,  where  "a"  Is  the  radius  of  the  cylinder. 

The  normalized  scattering  cross  sections  of  square  cylinders  are  shown 
In  Figures  4-23  through  4-27.  All  computed  results  are  normalized  by  irb, 
where  "b"  Is  one-half  the  width  of  the  square  cylinder  under  consideration. 
Twenty  expansions  have  been  used  in  all  computations  for  square  cylinders  of 
different  material  constants. 

Figures  4-28  through  4-30  show  the  characteristic  currents  (or  mode 
currents)  for  circular  cylinders  of  different  material  constants.  Fifteen 
expansion  functions  have  been  used  for  the  computation  of  mode  currents. 

For  representative  computations,  consider  a circular  cylinder  with 
ka  ■ 0.7  (where  "a"  is  the  radius  of  the  cylinder,  ct  ■ 9.5,  p{  ■ 1.0.  The 
contour  is  approximated  by  32  straight  lines  segments  of  equal  length  (the 
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line  segments  can  be  of  different  length),  and  15  expansion  functions  are  used 
for  both  electric  and  magnetic  surface  currents.  Figure  4-28  shows  the  char- 
acteristic currents  plotted  vs.  the  contour  length  variable  in  terms  of  a 
sequence  of  triangle  functions.  All  the  mode  currents  are  composite  currents. 
The  first  15  points  represent  the  electric  mode  current,  and  the  second  the 
magnetic  current. 

Figure  4-29  shows  the  characteristic  currents  for  a circular  cylinder 
with  er  * 50.0,  ur  - 1.0,  and  ka  ■ 0.7.  Figure  4-30  shows  the  characteristic 
currents  for  a circular  cylinder  of  er  - 2.56.  Note  that  every  mode  current 
is  normalized  by  its  maximum  magnitude. 

For  perpendicular  polarisation  (TE) , the  modal  solution  for  the  scattered 
field  agrees  extremely  well  with  the  scattered  field  computed  directly  from 
matrix  inversion.  The  scattering  cross  sections  using  characteristic  modes  are 
almost  identical  to  the  matrix  inversion  solutions  (the  differences  are  less 
than  0.001  db). 

To  be  specific,  Fig.  4-1  shows  the  normalized  scattering  cross  section  of 
a circular  cylinder  with  ■ 9.5,  ur  ■ 1.0,  and  ka  ■ 0.7  for  perpendicular 
polarization  (TE).  The  computed  scattering  cross  section  is  in  good  agreement 
with  harmonic  solution  (7].  The  maximum  deviation  is  0.65  db.  Figure  4-2 
glveo  the  normalized  scattering  cross  section  of  a circular  cylinder  with 
■ 20.0,  - 1.0,  and  ka  • 0.7,  for  perpendicular  polarization  (TE).  The 

maximum  deviation  from  exact  harmonic  series  solution  is  0.076  db.  Figure 
4-3  shows  the  normalized  scattering  cross  section  of  a circular  cylinder  with 
cr  ■ 50.0,  wr  ■ 1.0,  and  ka  • 0.7,  for  perpendicular  polarization  (TE).  Maximum 


/\ 
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deviation  from  exact  harmonic  solution  is  0.485  db.  The  scattering  cross 

section  shown  in  Fig.  4-4  is  for  a circular  cylinder  with  er  - 100.0, 

U ■ 0.01,  and  ka  ■ 0.7  for  perpendicular  polarization.  The  computed  solu- 

tion  is  in  excellent  agreement  with  the  exact  solution.  Maximum  deviation 

is  0.01  db.  Figure  4-5  gives  the  normalized  scattering  cross  section  of  a 

circular  cylinder  with  e ■ 1000.0,  w ■ 0.001,  and  ka  ■ 0.7,  for  perpen- 

r it 

dlcular  polarization.  Mote  that  the  computed  result  is  in  excellent  agree- 
ment with  the  calculations  of  a conducting  cylinder.  Maximum  deviation  is 
0.01  db.  The  conducting  cylinder  problem  can  be  viewed  as  a specialization 
of  the  more  general  material  cylinder  problem.  This  is  expected  to  be  true 
even  for  three-dimensional  objects.  Figure  4-6  shows  the  normalized  scatter- 
ing cross  section  of  a circular  cylinder  with  er  - 9.0,  nr  ■ 1.0,  and  ka  - 2.0, 
for  perpendicular  polarization  (TE).  Maximum  deviation  from  exact  harmonic 
solution  is  1.79  db.  Better  agreement  can  be  reached,  if  more  expansion  func- 
tions are  used.  The  scattering  cross  section  given  in  Fig.  4-7  is  for  a 
circular  cylinder  with  cf  » 9.0,  ur  ■ 1.0,  and  ka  ■ 1.0,  for  perpendicular 
polarization.  The  agreement  with  exact  solution  is  excellent.  Maximum  devi- 
ation is  0.013  db.  Figure  4-8  shows  the  normalized  scattering  cross  section 
of  a circular  cylinder  with  ■ 9.0,  uy  ■ 100.0,  and  ka  * 0.7,  for  perpen- 
dicular polarization.  Agreement  with  exact  solution  is  very  good.  Maximum 
deviation  ie  0.01  db.  Figure  4-9  shows  the  normalized  scattering  cross  sec- 
tion of  a circular  cylinder  with  er  ■ 9.0,  ■ 5.0,  and  ka  ■ 0.7,  for  per- 

pendicular polarization.  The  computed  result  is  in  good  agreement  with  exact 
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solution.  Maximum  deviation  is  0.3  db.  The  computed  scattering  cross  section 
of  a circular  cylinder  with  er  ■ 0.001,  » 1000.0,  and  ka  ■ 0.7,  for  per- 

pendicular polarization  is  shown  in  Fig.  4-10.  Note  that  the  cylinder  is  highly 
magnetic.  Maximum  deviation  from  exact  harmonic  solution  is  0.001  db.  The 
agreement  is  excellent.  Figure  4-11  shows  the  normalized  scattering  cross  section 
of  a circular  cylinder  with  * 1.0,  ur  ■ 1000.0,  and  ka  - 0.7,  for  perpendicular 
polarization.  Maximum  deviation  from  exact  solution  is  0,04  db.  Figure  4-12 
represents  the  computed  scattering  cross  section  of  a circular  cylinder  with 

e ■ 1.0,  u ■ 10.0,  and  ka  * 0.7,  for  perpendicular  polarization.  Maximum 
x r 

deviation  from  exact  solution  is  0.04  db.  Figure  4-13  shows  the  computed  scat- 
tering cross  section  of  a circular  cylinder  with  ■ 1.0,  ur  ■ 300,  and 
ka  - 0,7,  for  perpendicular  polarization.  Maximum  deviation  from  exact  hurmonic 
solution  is  0.2  db.  Figure  4-14  gives  the  normalized  scattering  cross  section 
of  a circular  cylinder  with  er  ■ 2.56,  ur  » 1.0,  and  ka  - 0.7,  for  perpendicular 
polarization.  Maximum  deviation  is  0.6  db.  Figure  4-15  Bhows  the  computed 

scattering  cross  section  of  a circular  cylinder  with  e ■ 1000.0,  u - 0.001, 

r r 

and  ka  « 0.7,  fox  perpendicular  polarization.  The  computed  solution  is  in 

excellent  agreement  with  exact  solution.  Maximum  deviation  is  0.01  db.  The 

computed  scattering  cross  sections  of  a circular  cylinder  with  ka  * 0.7,  arc. 

given  in  Fig.  4-16  for  three  different  sets  of  material  constants;  i)  * 1000.0, 

U ■ 1.0  1)  e ■ 10000.0,  u “ 1.0  iii)  e - 5.0,  “ 10~6.  All  are  for 

r r r r t 

perpendicular  polarization.  Figure  4-17  shows  the  normalized  scattering  cross 
section  of  a circular  cylinder  with  ■ 1000.0,  ut  ■ 0.001,  and  ka  ■ 0,7, 
for  parallel  polarization  (TM) . The  solution  agrees  excellently  with  conduct- 
ing cylinder  solution.  Maximum  deviation  is  0.023  db.  The  normalized  scattering 
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cross  section  of  a circular  cylinder  with  - 2.56,  ™ 1.0,  and  ka  - 0.7, 

for  parallel  polarization  is  3hown  in  Fig.  4-18.  Maximum  deviation  from  exact 
solution  is  05  db . Figure  4-19  represents  the  computed  scattering  cross  sec- 
tion of  a circular  cylinder  with  t ■ 20.0,  pf  ■ 1.0,  and  ka  ■ 0.7,  for  parallel 
polarization.  Maximum  deviation  from  exact  solution  is  0.2  db.  Figure  4-20 
shows  the  computed  scattering  cross  section  of  a circular  cylinder  with  » 50.0, 
V>T  *“  1>0,  and  ka  - 0.7,  for  parallel  polarization.  The  computed  solution  is  in 
excellent  agreement  with  exact  harmonic  solution.  Maximum  deviation  is  0.05  db. 
Figure  4-21  shows  the  computed  scattering  cross  section  of  a circular  cylinder 
with  er  ■ 4.0,  p ■ 1.0,  and  ka  ■ 0.7,  for  parallel  polarization.  Maximum  devi- 
ation from  exact  solution  is  0.2  db.  Figure  4-22  showB  the  computed  scattering 
cross  section  of  a circular  cylinder  with  er  * 9.5,  pf  ■ 1.0,  and  ka  * 0.7,  for 
parallel  polarization.  The  computed  solution  is  in  excellent  agreement  with 
exact  harmonic  solution  maximum  deviatwon  is  0.01  db.  Figure  4-23  shows  the 
computed  scattering  cross  sections  of  a square  cylinder  with  kb  ■»  1 .4,  for 
two  sets  of  material  constants:  i)  er  » 1000.0,  p^  - 0.001  li)  er  ■ 1000.0, 

pf  ■ 1.0,  all  for  perpendicular  polarization.  For  square  cylinders,  there 
are  no  exact  solutions.  Figure  4-24  shows  the  computed  scattering  cross 
section  of  a square  cylinder  with  er  ■ 10000.0,  pf  * 0.0001,  and  kb  ■ 1.4, 
for  perpendicular  polarization.  The  computed  solution  has  been  compared  with 
the  solution  of  a conducting  square  cylinder  by  using  E-field  formulation  [13]. 
Maximum  deviation  is  0.1  db.  Figure  4-25  showB  the  computed  scattering  cross 
section  of  a square  cylinder  with  ■ 9.0,  pr  • 1.0,  and  kb  ■ 1.4,  for  per- 

pendicular polarization.  Figure  4-26  shows  the  computed  scattering  cross 
section  of  a square  cylinder  with  ■ 100.0,  p^  » 1.0,  and  kb  - 1.4,  for 
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parallel  polarization.  Figure  4-27  shows  the  scattering  cross  section  o£  a 
square  cylinder  with  er  - 10000.0,  y^  - 0.0001,  and  kb  «■  1.4,  for  parallel 
polarization.  The  computed  result  is  in  excellent  agreement  with  conducting 
square  cylinder  solution.  Figure  4-28  shows  the  lowest  order  characteristic 
currents,  plotted  as  a function  of  the  contour  variable.  The  currents  are 
normalized  by  choosing  their  maximum  amplitude  to  be  unity.  The  characteristic 
currents  are  for  a circular  cylinder  with  m 9.5,  y^  - 1.0,  and  ka  ■ 0.7,  for 
perpendicular  polarization.  The  electric  part  of  each  characteristic  current  is 
circumferentially  directed  and  the  magnetic  part  is  axially  directed.  The 
scattering  cross  section  computed  from  modal  solution  is  almost  identical  to 
that  from  matrix  inversion.  Figure  4-29  shows  the  normalized  characteristic 
currents  for  a circular  cylinder  with  ■ 50.0,  yf  ■ 1.0,  and  ka  ■ 0.7,  for 
perpendicular  polarization.  Figure  4-30  gives  the  normalized  characteristic 
currents  for  a circular  cylinder  with  er  « 2,56,  y^  » 1.0,  and  ka  ■ 0.7,  for 
perpendicular  polarization. 

The  purpose  of  this  work  is  to  show  the  feasibility  that  a surface 
formulation  for  the  theory  characteristic  modes  can  be  applied  to  the  solu- 
tion of  scattering  from  material  objects.  For  large  cylinders,  more  expan- 
sion functions  are  needed.  No  attempt  has  been  made  to  treat  large  objects. 

It  is  expected  that  this  is  one  of  the  important  areas  for  future  research. 

Many  questions  are  still  left  unanswered  in  the  interpretation  and  applica- 
tion of  characteristic  modes  to  material  objects.  It  is  hoped  that  this  work 
will  be  of  some  value  to  future  researchers  in  their  effort  to  gain  a complete 
understanding  of  the  theory  of  characteristic  modes. 
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The  eigenvalue  equation  (3-61)  is 

[T2][f]  - X [T^ff] 


(4-1) 


and  the  expression  for  the  Rayleigh  quotient  associated  with  equation 


(4-1)  is 


[£iltT21Ifi:i 

xi  " IfilTii'lTfJ 


(4-2) 


The  computed  eigenvalues  and  their  corresponding  eigenvectors  should  satisfy 
equation  (4-2).  The  Rayleigh  quotient  check  is  important  because  it  gives 
some  -verification  to  the  approximations  used  in  numerical  computation. 

The  quadratic  term  [f^] [T^] [f^)  deserves  some  elaboration  since  it 
qppears  frequently  in  equations.  Note  that 

IfHTjHf]  - hOtunTjtuHx] 

- uitunx] 

- (x2l  (4-3) 

It  has  already  been  pointed  out  in  Chapter  3 that  [x2]  is  the  component  of 
an  eigencurrent  f that  lies  within  the  null  space  of  T^,  in  other  words, 

[x2]  does  not  radiate.  Since  approximations  are  made  in  the  computational 
procedures,  the  elgencurrents  will  not  be  absolutely  exact.  Consequently, 
the  second  quadratic  term  on  the  right  hand  side  of  equation  (4-3)  will 
differ  from  zero,  but  it  should  be  much  smaller  than  the  first  quadratic 
term.  To  a certain  degree,  this  will  give  some  indication  of  the  accuracy 
of  the  computed  elgencurrents.  The  first  quadratic  term  at  the  right  hand 
side  of  equation  (4-3)  can  be  further  expressed  as 
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[Xj]  [x^ 

- tw^1/2]  [y]  lu^  [y] 

- [Tlh^H^Hp^Hy] 

- [y][y] 

-I  (if  {y^>  are  orthonormal)  (4-4) 

In  numerical  computation,  approximations  are  inevitable.  Some  special 
analytical  manipulations  such  as  those  discussed  above  can  often  provide 
added  insight  to  the  correctness  of  the  numerical  results. 
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Exact 

This  method,  20  expansion  functions. 
Max.  dev.  0.01  db 


$ (degrees) 

Jig.  4-9,  Noraalized  scattering  crosB  section  of  a circular  cylinder  with 
€ ■ 9.0,  J*  ■ 5.0,  ka  - 0.7,  perpendicular  polarization  (TE) . 
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Exar-  r 


*i Ills  method,  20  expansion  functions. 
Max.  dev.  0.04  db. 
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20  40  60  80  100  120  140  160  180 

4>(degrees) 

Fig.  4-12.  Normalized  scattering  cross  section  of  a circular  cylinder  with 
t'r  ■ 1.0,  - 10.0,  ka  ■ 0.7,  perpendicular  polarization  (TE) . 
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Exact 


This  method,  20  expansion  functions. 
Max.  dev.  0.2  db. 


I I I I I I I I 1 1 

20  40  60  80  100  120  140  160  180 

$(degrees) 

Fig.  4-19.  Normalized  scattering  cross  section  of  a circular  cylinder 

with  - 20.0,  « 1,0,  lea  - 0.7,  parallel  polarization  (TM) . 
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DISCUSSION 

A surface  formulation  la  developed  for  solving  two-dimensional 
electromagnetic  scattering  problems.  A basic  theory  for  characteristic 
modes  of  dielectric  and  magnetic  bodies  based  on  the  surface  formulation 
is  derived.  The  method  of  computing  characteristic  modes  can  be  used  for 
homogeneous  material  bodies  of  arbitrary  shape  provided  the  body  1b  not 
electrically  large.  The  characteristic  modes  of  material  bodies  have  most 
of  the  properties  of  those  for  perfectly  conducting  bodies,  and  should  find 
similar  uses.  The  theory  presented  here  is  In  contrast  to  that  for  the 
volume  formulation  [4].  The  basic  difference  is  that  the  current  in  the 
material  body  has  been  treated  as  equivalent  surface  currents  Instead  of  a 
volume  distribution.  The  characteristic  currants  are  real  and  their  cor- 
responding eigenvalues  are  also  real.  The  eigenvectors  given  by  equation 
(3-73)  are  those  corresponding  to  the  lowest  eigenvalues,  and  they  are 
usually  very  efficient  radiators.  Characteristic  currents  associated  with 
large  eigenvalues  generally  Indicate  higher  order  modes  which  do  not  radiate 
very  much. 

Two  ways  for  computing  the  scattered  fields  are  given  here.  The  simple 
material  cylinders.  The  matrix  inversion  method  is  easier  to  use  and  gives 
very  good  results.  The  characteristic  mode  method  may  require  slightly 
longer  computing  time,  but  it  does  provide  more  insight  into  the  problem. 

As  in  the  conducting  body  case,  the  characteristic  mode  method  should  prove 
to  be  of  value,  both  theoretically  and  computationally  for  scattering  and 
radiation  problems.  The  versatility  of  characteristic  modes  has  been 
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adequately  demonstrated  in  analysis  and  synthesis  problems  dealing  with 
conducting  bodies.  The  two  approaches  are  based  on  a surface  formulation, 
and  they  require  the  material  body  to  be  homogeneous  since  the  unknowns 
are  surface  currents.  For  Inhomogeneous  bodies  the  surface  formulation  is 
not  appropriate,  and  a volume  current  distribution  must  be  used  which  re- 
quires sample  points  Inside  the  scattering  body. 
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APPENDIX  A 

MATRIX  ELEMENTS  FOR  PARALLEL  POLARIZATION 


For  the  incident  field 


E ■ u e~4fe  (A-l) 

***•  z 

VA 

the  following  formulas  are  obtained  ( u i«  the  axially  directed  unit  vector) . 

E 

The  procedures  involved  are  identical  to  that  given  in  Chapter  2,  except  that 
the  directions  of  the  surface  currents  are  different. 

A. 1 Formulas  for  [21  Matrix  Elements 


where 
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Z ■ At  At  H^(kR  ) (non-coincident  intervals) 

p q o pq 

2 YkAt 
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A. 2 Formulas  for  [Bl  Matrix  Elements 
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A. 3 Formulas  for  [C]  Matrix  Elements 
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A.4  Formulas  for  [Y]  Matrix  Elements 
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A.5  Excitation  Matrix  Elements 
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l I I.MNf  INUf 

**f  *»1*H 
l*.  r,n*it|Nuf 
■ I *0 

1 Id*. i 

IM.MRI.IM  IH.R.IH 
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iisn«v*((i.nni  \ if>m~v*n,nnn?nfm  m 1 1 1 
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